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AI Bao-Quan, WANG Xian-Ju, LIU Guo-Tao and LIU Liang-Gang * 
Department of Physics, ZhongShan University, GuangZhou, P. R. China. 

Nakano M 

Department of information science. University of Occupational and Environmental 

Health, Japan. 
Matsuura H 

Department of Project Center, National Graduate Institute for Policy Studies, Japan. 

The SIRM (Stochastic Inchned Rods Model) proposed by H. Matsuura 
and M. Nakano can explain the muscle's motion perfectly, but the inter- 
molecular potential between myosin head and G-actin is too simple and 
only repulsive potential is considered. In this paper we study the SIRM 
with different complex potential and discuss the effect of the spring on the 
system. The calculation results show that the spring, the effective radius of 
the G-actin and the intermolecular potential play key roles in the motion. 
The sliding speed is about 4.7 x W~^m/s calculated from the model which 
well agrees with the experimental data. 
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I. INTRODUCTION 

The great progress have been made on the observations of molecular motor since 
H. F. Huxley and H. E. Huxley [1] [2] [3] proposed the rotating cross-bridge model to 
study muscular mobility in 1954. More and more biologists and theoretical biophysicists 
concentrate on studying the mechanism of the living system and mincing the biological 
structure to invent subtle artificial instrument. Based on different system, three different 
families of motor proteins have been identified [4] [5] [6], Kinesins and dyneins move along 
tubulin filaments [7], myosin move along actin filaments [8] [9]. The motion mechanism of 
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the three kinds of motors can be described as follows [4] [10] [11]: firstly, the motor 
binds Adenosinetriphosphate (ATP), then hydrolyzes the bound ATP and absorbs the 
energy, subsequently, the motor releases the products Adenosinediphosphate (ADP) and 
Phosphate (P), so the motor move constantly within the chemical cycle. The molecular 
motors play a key role in transuding chemical energy into mechanical work at a molecular 
scale. Up to date there have been three kinds of archetype models proposed to study the 
molecular motor system: namely the fluctuation models, the phase transition models and 
the resonance models. 

One archetype fluctuation model is ratchets model. This model looks molecular motors 
as heavily damped Brownian particles and the model explains the molecular motor's 
motion by studying a single particle in a potential inspired by Feynman's ratchet. M. O. 
Magnasco [12] and R. D. Astumian [13] [25] have proposed the forced thermal ratchets 
model, the fluctuant force break the symmetry of the system, so that the system can 
absorb the energy from the surroundings, the motor undergoes similar Brownian motion 
and produces net transport. Although the fluctuation models have discussed different 
aspects of molecular motors, there still exits a matter of debate on the mechanism of its 
motion. 

The second typical model is phase transition model: the two state model proposed by 
F.Julicher [4] and isothermal ratchet model proposed by A. Pameggiani [14]. Using the 
basic idea of the Ising model for phase transitions, they have studied a single motor or 
many motors moving in asymmetry periodic potential. In these models the microtubule 
and molecular motor are treated as periodic polar railway tracks and carnot engines with 
cargoes, the motor consumes chemical energy and causes the system conformation change, 
so the phase transition occurs between two states. The models give a good explanation 
on the efficiency of the motor. 

The newly type resonance model is named the Stochastic Inclined Rods Model (SIRM) 
proposed by H. Matsuura and H. Matsuura [?], in the model the energy of motion was 
supplied from the random noise and the system always moves to one direction by using 
stochastic resonance. The movement of the system does not break the second law of 
thermodynamics, because the actin-myosin system is open to the surroundings and the 
energy flows in from their surroundings. The SIRM presents a perfectly sliding mechanism 
of the actin-myosin system, but the interaction is too simple and only repulsive force is 
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considered. The main aim of our paper is tliat based on SIRM we study muscle's motion 
under a complex intermolecular potential and discuss the effect of the spring on the 
system. 



II. MODEL AND FORMALISM 



In order to describe the motion of actin-myosin system, we simulate it by a mechan- 
ical model as shown in Fig. 1. SIRM is composed of three parts: an inclined spring, a 
myosin head and myosin bundle. Our model proposed a microscopic mechanism for the 
actin-myosin system: The random noise is from the water molecules generated through 
the heat energy of ATP hydrolysis. ATP supphes the heat energy and provides heat 
energy to the surrounding water molecules, which accelerates the movement of the water 
molecules. So the random noise (the water molecules) can interact with the myosin head, 
and the resonance occurs between the myosin head and the noise, then it collides with 
a G-actin and obliquely kicks the G-actin sphere, because the direction of vibration is 
inclined against the line of the actin fibers, myosin molecules obtain the propellant force 
along the direction of the fiber and the myosin can move to one direction. 



Let us set up the dynamic equations of the system. As for the intermolecular potential 
between myosin head and G-actin we adopt Jennard- Jones potential: 



Where ri — y{x — XiY + {y — UiY — R, {xi, yi) means a center of i-th G-actin glove, 
and R is its effective radius, p, q, Ul^,jk are the parameters of the potential. 
The potential of the myosin rod is approximately expressed as follows: 



Fig. 1 
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Where Kq, Ki are tangent constant and radial constant of the spring, and L is initial 
length of the spring , tan(^) = , 6', 6*0 are current angle and initial angle between the 

rod and level axis. (a;2, ^2)18 the center of gravity of myosin bundle. 

The equations of the myosin head are: 

d'y d{Ua + Us) , ^dy 

With regard to the center of gravity {x2, ^2) of the filament, we set a similar equation: 

(9^X2 dUs , . dx2 , . 

The variable t/2 is fixed since the myosin filament does not significantly move along the 
^/-direction compared to the x-direction. 



Where a, (3, rj are viscous constants, M is the mass of the myosin bundle and m is 
the mass of the myosin head. Fx{t), Fy{t), F^^it) are fluctuation of the thermal noise, the 
noise is Gauss white noise and it satisfies the following fiuctuation-dissipative relations 
[21] [22] [23]. 

< >= 0. (6) 

< F,(t)F,(s) >= 2kBTC5a,b5(t - s). (7) 

where a,b — x,y,X2;C — 0^,(3,7]; ks is Boltzmann constant. T is absolute temperature, 
t, s are time. 

III. CALCULATION RESULTS 

In order to solve the above equations we adopt a numerical method. The parameters 
of the equations are shown in the table. 
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Table. 1. The Parameters of the Equations. 



1 unit time 


10-^s 


a,(3,r] 


1 




15000 


1 unit length 


lO^^m 


p,q 


100 


m 


80 


1 unit mass 


1.66 X 10-2^kg 


L 


10 


M 


400 


Ke 


450000 


Ki 


4500 


Oo 


45° 



We solve the equations (l)-( 7) with a numerical method and the numerical results are 
shown in Fig.2-Fig.6. 



Fig. 2 

The dashed line shows the position of myosin filament and the solid line describes the 
displacement of myosin head. Prom the figure we can know that the head and the fil- 
ament, bound by the spring, move forward along the same direction as a whole. The 
movement of the head is a translational one as a whole, though the trace has irregularity 
and randomness like Brownian particles. Wc can know from the figure the average gliding 
speed of the myosin bundle is 4.7 x 10~^m/s, On the other hand, the experimental data 
by Yanagida et al. is that the maximum gliding speed is 7 x 10~^m/s. It is impressive 
that our prediction well agree with the experimental data. 



Fig. 3 

The solid line shows the stretching motion of the myosin rod, the vertical displacement 
of myosin head is described with the dashed line. The figure shows that the collision 
between the G-actin and the myosin head distorts their relative motion from a trigono- 
metric function. This trace of the motion is given by an accurate trigonometric function 
if there are no interactions or the collisions. Truly the strain of the trigonometric function 
indicates that the myosin head kicks obliquely the G-actin when the spring is stretched 
by the resonance vibration. 



Fig. 4 



5 



From the figure we can know that the constant of the spring play a key role in the mo- 
tion, if the constant is zero, the system can not move. And if the constant increases, the 
velocity of the SIRM will increase rapidly. 



Fig. 5 

The figure can show that the effective radius of the G-actin is important for the system 
to move. If the G-actin is flat {R — 0) the SIRM cannot move, on the other hand, if the 
effective radius of the G-actin is too big {R — Anm) the system cannot move, either. The 
system can move at the highest speed R — 2.bnm. 



Fig. 6 

From the flgure we can know that the repulsive potential can accelerate the system's mo- 
tion and the potential between the G-actin and myosin can decide the motion's velocity. 



IV. SUMMARY AND CONCLUSION 



Based on the SIRM we study the movement of actin-myosin system with a Jennard- 
Jones potential. From the numerical results we can know there is relative sliding motion 
between the filament and the G-actin and the average gliding speed is 4.7 x 10~^m/s 
which well agree with the experimental data. When the system absorbs the energy from 
the thermal noise constantly through stochastic resonance, the intermolecular potential 
and the inclined rod make the filament move to one direction. The system can convert 
the random noise to one directional motion. The effective radius of G-actin, the constant 
of the spring and the intermolecular potential are important for the system to move. If 
the constant of the spring is small, the system cannot move, and if the effective radius of 
G-actin is too small or too big the system cannot move, either. The repulsive potential 
makes the system move fast. From the results we can predict the state of the muscle: 
if the constant of the spring or the effective radius of G-actin become small, the SIRM 
can not move at all and muscle is fatigue, if the effective radius of G-actin is too big, the 
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velocity of the system is zero and the actin-myosin system may be breakdown. 

SIRM is thermally open to the outer surroundings and it has a outer heat source of 
ATP. The heat or energy consumed by SIRM is directly supplied from the surroundings 
(i.e. the thermal motions of water molecules). The entropy always increases since the 
flow is positive, therefore, the system does not break the second law of thermodynamics. 
There have been many discussions on the movement of ratchets in a heat bath. Among 
them, we only consider one famous discussion of Feynman [24]. Feynman showed from 
statistical discussions that the ratchets under an isothermal temperature cannot turn as a 
total since the ratchets turn to both sides with equal probability. However, he also showed 
in the same section, 46-2, of Ref. [24] that it can turn to one side if a load is suspended 
from the one side even in an isothermal bath. Our model can work even under isothermal 
temperature conditions, since it has a load, which is an intermolecular potential of the 
G-actin fastened to a Z-membrane. 

Noise in dynamical systems is usually considered a nuisance. But in certain nonlinear 
systems including electronic circuits and biological sensory apparatus the presence of noise 
can in fact enhance the detection of the weak signals. From the above results, we point 
out the possibility that the actin-myosin system can use thermal noise around for the 
movement in similar way. This viewpoint is completely new as the traditional viewpoint 
thinks the water disturbs the motion as viscosity. 
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V. FIGURE AND TABLE CAPTIONS 



Table. 1. The parameters of the equations, a, (3,7] arc viscous constants, M,m are 
the mass of myosin bundle and myosin head , respectively, L is initial length of the spring, 
p,q,U^,jk are the parameters of the potential, Kq^Ki are tangent constant and radial 
constant of the spring. 

Fig. l.The model for actin- myosin system. The coordinates of the head of myosin 
and the end of myosin connected to the filament are (x, and {x2,y2), respectively. The 
body of the myosin is assumed to be represented as a spring. 

Fig. 2. The dashed hne shows the position of myosin bundle and the sohd hne describes 
the displacement of the myosin head. The parameters are shown in the figures also. 
Fig. 3. The dash line shows the vertical movement of the myosin head and the solid line 
presents the stretching displacement s of the rod and s = \J [x — x^i)^ + (l/ — — -f'- 
Fig. 4. The relation between the constant of the spring and the x-displacement of the 
myosin bundle is shown in the figure. Prom the up to the bottom, different lives corre- 
spond to Ki: 5000,3000,1000,0, respectively. 

Fig. 5. The relation between the effective radius of the G-actin and the x-displacement 
of the myosin bundle is shown in the figure. From the up to the bottom, different lives 
correspond to i?: 2.5, 1, 4, 0, respectively. 

Fig. 6. The figure indicates the relationship of the intermolecular potential and the x- 
displacement of the myosin bundle. Prom the up to the bottom, different lives correspond 

to ik: 4, 3, 2, 1, respectively. 
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